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Abstract 

This paper is concerned with the time-step condition of commonly-used linearized 
semi-implicit schemes for nonlinear parabolic PDEs with Galerkin finite element approx- 
imations. In particular, we study the time-dependent nonlinear Joule heating equations. 
We present optimal error estimates of the semi-implicit Euler scheme in both the L 2 
norm and the H 1 norm without any time-step restriction. Theoretical analysis is based 
on a new splitting of error function and precise analysis of a corresponding time-discrete 
system. The method used in this paper is applicable for more general nonlinear parabolic 
systems and many other linearized (semi)-implicit time discretizations for which previ- 
ous works often require certain restriction on the time-step size t. 
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1 Introduction 

In the last several decades, numerous effort has been devoted to the development of efficient 
numerical schemes for nonlinear parabolic PDEs arising from a variety of physical appli- 
cations. A key issue to those schemes is the time-step condition. Usually, fully implicit 
schemes are unconditionally stable. However, at each time step, one has to solve a system 
of nonlinear equations. An explicit scheme is much easy in computation. But it suffers 
the severely restricted time-step size for convergence. A popular and widely-used approach 
is a linearized (semi)-implicit scheme, such as linearized semi-implicit Euler scheme. At 
each time step, the scheme only requires the solution of a linear system. To study the 
error estimate of linearized (semi)-implicit schemes, the boundedness of numerical solution 
(or error function) in L°° norm or a stronger norm is often required. If a priori estimate 
for numerical solution in such a norm cannot be obtained, one may employ the induction 
method with inverse inequality to bound the numerical solution, such as 

- t/^IUoo < CJ/t-^II^C-,^) - Z7^|| i2 < C^-^ 2 ^ + (1.1) 
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where u(-,t n ) and Ufi are the exact solution and numerical solution, respectively, Rh is 
some projection operator and d is the dimension. The above inequality, however, results 
in a time-step restriction, particularly for problems in three spatial dimensions. Such a 
technique has been widely used in the error analysis for many different nonlinear parabolic 
PDEs, e.g., see [H H HBJ EQl EI] for Navier-Stokes equations, [21 HUES] for nonlinear Joule 
heating problems, [151 [23 [2Z] for porous media flows, [3 [T2J CGE [2H] for viscoelastic fluid 
flow, [221 [35] for KdV equations and [TO], [29] for some other equations. In all these works, 
error estimates were established under certain time-step restrictions. We believe that these 
time-step restrictions may not be necessary in most cases. In this paper, we only focus 
our attention to a time-dependent and nonlinear Joule heating system by a linearized semi- 
implicit scheme. However, our approach is applicable for more general nonlinear parabolic 
PDEs and many other time discretizations to obtain optimal error estimates unconditionally. 
The time-dependent nonlinear Joule heating system is defined by 

^-Au = a(u)\V ( j ) \ 2 , (1.2) 

- V • (<t(u)V0) = 0, (1.3) 

for x € f2 and t G [0, T], where Q. is a bounded smooth domain in M. d , d = 2,3. The initial 
and boundary conditions are given by 

u(x,t) = 0, 4>{x,t) = g(x,t) for x £ d£l, t £ [0,T], 
u(x,0) = uq{x) for x £ £1. 

The nonlinear system above describes the model of electric heating of a conducting body, 
where u is the temperature, (f> is the electric potential, and a is the temperature-dependent 
electric conductivity. Following the previous works [11[ 136]. we assume that a £ W ,0 °(R) 
and 

k < a(s) < K (1.5) 

for some positive constants k and K. 

Theoretical analysis for the Joule heating system was done by several authors [3 
EH EH E21 E3]. Among these works, Yuan [33] proved existence and uniqueness of a C a 
solution in three-dimensional space. Based on this result, further regularity can be derived 
with suitable assumption on the initial and boundary conditions. Numerical methods and 
analysis for the Joule heating system can be found in [21 [U HH EQl EU E7J EH] . For the system 
in two-dimensional space, optimal I? error estimate of a mixed finite element method with 
the linearized semi-implicit Euler scheme was obtained in [36] under a weak time-step 
condition. Error analysis for the three-dimensional model was given in [11], in which the 
linearized semi-implicit Euler scheme with a linear Galerkin FEM was used. An optimal 
L 2 -error estimate was presented under the time step restriction r < 0{h 1 / 2 ). A more 
general time discretization with higher-order finite element approximations was studied in 
[2]. An optimal L 2 -norm error estimate was given under the conditions r < 0(h 3 / 2p ) and 
r > 2 where p is the order of the discrete scheme in time direction and r is the degree of 
piecewise polynomial approximations used. No optimal error estimates in i? 1 -norm have 
been obtained. 

The main idea of this paper is a splitting of the numerical error into the temporal 
direction and the spatial direction by introducing a corresponding time-discrete parabolic 
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system (or elliptic system). Error bounds of the Galerkin finite element methods for the 
time-discrete parabolic equations in certain norm is dependent only upon the spatial mesh 
size h and independent of the time-step size r. If a suitable regularity of the solution of 
the time-discrete equations can be proved, numerical solution in the L°° norm (or stronger 
norm) is bounded unconditionally by the induction assumption 

\\R h U n - U r ^\\ L oo < Ch- d ' 2 \\R h U n - U%\\ L 2 < Ch- d ' 2 h r+1 (1.6) 

where U n is the solution of the time-discrete equations. With the boundedness, optimal 
error estimates can be established for the fully discrete scheme without any time-step restric- 
tion. In this paper, we analyze the linearized (semi-implicit) backward Euler scheme with 
the standard Galerkin approximation in spatial directions for the nonlinear Joule heating 
system (|1 .2j) - f)l .4[) . With the splitting, we present unconditionally optimal error estimates 
in both the L 2 norm and the H 1 norm. 

The rest of the paper is organized as follows. In Section 2, we present the linearized 
semi-implicit Euler scheme with a linear Galerkin finite element approximation in spatial 
direction and our main results. After introducing the corresponding time-discrete parabolic 
system, we provide in Section 3 a priori estimates and optimal error estimates for the time- 
discrete solution, which imply the suitable regularity of the time-discrete solution. With 
the regularity obtained, we present optimal error estimates of the Galerkin finite element 
solution in L 2 -norm without any time-step restriction, and the optimal error estimate in 
H 1 norm follows immediately due to the nature of our approach. The concluding remarks 
are presented in Section 4. Extension to r-order Galerkin finite element approximation is 
straightforward with the corresponding assumptions of regularity. 



2 Galerkin methods and main results 

Let Q, be a bounded convex and smooth domain in M. d (d = 2,3). For any integer k > 
and 1 < p < oo. Let Wr(f2) be the Sobolev space with the norm 



where 



\P\<k 



dx^ 1 ■ ■ ■ dx d d 

for the multi-index (3 = (/3 1; • • • , (3d), f3\ > 0, • • • , fid > 0, and \/3\ = fi\ + ■ ■ ■ + /3d- For any 
integer k > and < a < 1, let C k+a (Q) denote the usual Holder space with the norm 

ii/ii^-Eii^/ii cm + E^p ™^ 

\/3\<k \/3\=k x ' yen 1 yl 

and let Cb(f2) be the space of continuous functions on f2 vanishing on the boundary 90. 
For any Banach space X and function / : [0, T] — > X, we define the norm 



T 

LP((0,T);X) 



n n/(*)ii^y, i<p<oo, 

esssup t6 ( 0)T )||/(t)||x, p = oo. 
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With the boundary conditions (jl.4p . the weak formulation of the system (|1.2|) - (|1.3|) is 
defined by 

(u t , £„) + (Vn, V£„) = (<7(u)|V^ 2 , (2.1) 
(<t(«)V0, V^) = (2.2) 

for any £ u ,£<^ G -Ho(^) an( ^ a,e - * e (0>^)- 

Let TTh be a regular division of f2 into triangles in M 2 or tetrahedras in M 3 , i.e. Q = Lljttj, 
and denote by h = maxjjdiam ttj} the mesh size. For a triangle Qj at the boundary, we 
define Qj as the triangle with one curved side (or a tetrahedra with one curved face in ]R 3 ) 
with the same vertices as Qj, and set Dj = For an interior triangle, we set Qj = Qj 

and Dj = 0. For a given division 7r/j, we define the finite element spaces [26J: 

Vh = {vh G C(O) : Vhlfij is linear for each element and v h = on -Dj}, 
Sh = £ C(f2) : vh\q. is linear for each element}. 

It follows that Vh is a subspace of Hq(£1) and is a subspace of H 1 ^). For any function 
v G S/j, we define AhV as the function which satisfies AhV = on Dj and AhV = v on Tj. 
We define 11^ : C(Q) — > Sh to be the Lagrangian interpolation operator, i.e. HhV coincides 
with v at each vertex of the triangular division of and set Tlh = AhUh- Clearly, Uh is a 
projection operator from Co(£l) onto Vh- 

Let {t n }n =0 be a partition in the time direction with t n = nr, T = Nt and 

u n = u(x,t n ), 4> n = (j)(x, t n ) . 

For any sequence of functions {/ n }^L 0) we define 

fn+l J"n 

D t f n+1 = - — . 

r 

For simplicity, we assume that g G H 1 ^). The fully discrete finite element scheme is 
to find U», $1 - g n G V h for n = 0, 1, • • • ,N such that for all £ U) ^ G F h 

(A^ n+1 , e«) + (Vt£ +1 , Ve«) = H^)|V^| 2 , Q, (2.3) 
(a(^)V^, V^) = 0, (2.4) 

with the initial conditions U? = RhU°. 

In the rest part of this paper, we always assume that the solution to the initial/boundary 
value problem (|1.2p - (|1.4[) exists and satisfies 

IMIl°°((0,T);H 2 ) + IMIl°°((0,T);L 2 ) + IK \\l 2 {(0,T);H 2 ) + \\ u tt ||l 2 ((0,T);L 2 ) + \\ u o\\h 2 
+ H0lli«>((O,T);W a . 12 /8) + ||^t||l,2((0,r);Hl) + II WIIl^O.T);^) 

+ II#IIl«=((0,T);W 2 . 12 / 5 ) + llfl't|lL a ((0,T);H 1 ) + \\Vg\\L°°((o,Ty,c°>) < C. (2.5) 

We denote by C a generic positive constant, which is independent of n, h and r and e a 
generic small positive constant. We present our main results in the following theorem. 

Theorem 2.1 Suppose that the system ( U.^j )-( TOj) i/ie initial and boundary conditions 
has a unique solution (u,<f>) satisfying \2. 5\) . Then there exist positive constants tq 



4 



and h,Q such that when t < tq and h < ho, the finite element system 12.3\) - [K4\ ) admits a 
unique solution (Ufi, $^), n = 1, • • • ,N, such that 

max \\UJi - u n \\ L2 + max ||^ - ^|| L2 < C(r + h 2 ), (2.6) 

l<n<yV l<n<N 

max \\U£ - u n \\ m + max ||d>™ - 4> n \\m < C(r + h). (2.7) 

l<n<N l<n<N 

For £/° = no and $°, we define C/ n and <I ) ' 1 to be the solution of the following discrete 
parabolic system (or elliptic system) 

D t U n+1 - AU n+l = o-(U n )\V$ n \ 2 , 0<n<N-l, (2.8) 

- V • ((i([/")V$ n ) = 0, < n < N, (2.9) 

with the boundary conditions 

U n+1 (x) = 0, $ n (x) = g{x,t n ) for x € 8CI. (2.10) 

We will present the proof of Theorem 12. II in the next two sections. The key to our proof 
is the following error splitting 

\\Uh - uTl \\ < ll e l + Kll + \\ un ~ Rh.U n \\, 

\M ~ < \\v n \\ + \\Vh\\ + 11$'" " Ph$ n \\ 
for any norm || • ||, where 

e n = U n - u n , e n h = XJl- R h U n , 

with Rh : Hq(Q) — > being the Riesz projection operator, i.e. 

(V(u - R h v), Vw) = 0, for all v G H^(Q.) and w G V h . 

and P£$ n = g(-,t n ) + U h ($ n - g{-,t n )) for n = 0, 1, 2, • ■ ■ , N. 

With the definition of the operator U^, Hh, P% and Rh, the following estimates hold 
[23]: for any 2 < p < oo, there exists a positive constant C (independent of the function v) 
such that 

\\v - Ilhv\\LP + h\\v - ILhvWw 1 - Ch 2 \\v\\ W 2, ( 2 -ll) 
\\v - U h v\\ LP + h\\v - Il h v\\ w i < Ch 2 \\v\\ W 2, (2.12) 
||$« _ p«$«|| LP + h\\$n _ P^ n \\ w i < Ch 2 \\$ n - g n \\ W 2, (2.13) 
\\v - R h v\\ LP + h\\v - R h v\\ w i < Ch 2 \\v\\ W 2 (2.14) 
hold for all v G W 2 (Q). 

3 Error estimates 

We analyze the error function (e n , rj n ) from the linearized semi-implicit Euler scheme (time- 
discrete system) and the errors function (e^,??^) of the Galerkin finite element method for 
the time-discrete system in the following two subsections, respectively. 
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3.1 The time-discrete solution 

In this subsection, we prove the existence and uniqueness of the time-discrete system 
(|2.10p and establish the error bounds for (e n ,rj n ). 

Theorem 3.1 Suppose that the system il.2\) - (L4\ ) has a unique solution (u,<f)) satisfying 
\2. 5\) . Then there exists a positive constant tq such that when r < tq, the time- discrete 
system \2.8\) - l2~W\) admits a unique solution (U n ,& n ) such that 



(3.1) 



/ N \ V 2 

max \\U n \\ H 2 + max \\D t U n \\ L 2 + ( V rll D t U n f m | < C, 

l<n<N l<n<N J 

max \\<$> n \\ H 2+ max \\V<5> n \\ L ™ < C (3.2) 

Kn<JV Kn<N 



and 



max ||e"|| H 'i+ max ||r? n ||rfi < Ct. C\'\\ 

l<n<N U UH l<n<N lU ~ \ 6 ' 6 ) 

Proof We rewrite the system ()1.2[) - ()1.3|) by 

D t u n+1 - Au n+1 = a(u n )\V(j) n \ 2 + B% +1 , (3.4) 

- V • (a(u n )V<j> n ) = 0, (3.5) 

where is the truncation errors due to the time discretization, i.e. 



du 



M* 1 = D t u n+L - ^ + {a{u n+l ) - o{u n ))\V(j) n+1 \ 2 

ut t—t n +i 

+ ( j(n n )V(0 n+1 + <f> n ) ■ V(0 n+1 - <p n ). 
With the regularity given in (I2.5p . we have 

< C, En=o \\Ri +1 \\hr < Cr\ (3.6) 
Subtracting the equations (13.4p - (|3.5p from the equations (12.8p - (|2.9p . respectively, we obtain 

D t e n+1 - Ae n+1 = (a(U n ) - o{u n ))\V<p n \ 2 

+ a{U n )(V(j) n + V$ n ) • Vrf + i^ +1 , (3.7) 

- V • (a(U n )Vr] n ) = V • [(a(u n ) - a{U n ))V^ n ] , (3.8) 
with the boundary condition e n+1 = n n = on dQ,. An alternative to the last equation is 

- V • {a(u n )Vif) = V • [{a{u n ) - <j{U n )){V<i? 1 + Vrf 1 )}. (3.9) 
Multiplying the equation (I3.8P by rf l+l and integrating the result over we have 

HVrfHi^CinMlVT/l^ 

which leads to 

||V771l 2 <C\\e n \\ L 2. (3.10) 
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Similarly, multiplying (|3.7p by e n+1 and integrating it over Q gives 
AQ||e" +1 ||i 2 )+||Ve^||i 2 

< C|Ie n || i a||e n+1 || La |lV0 n ||L«» + {a{U n ){V cj) n + V$ n )e n+1 , Vr/ n ) 

, II pn+l II „ Ho n + 1 II „ 
+ \\ K 1 \\L 2 \\ e IlL 2 • 

By (|2.9|) and using integrating by part, 

\{a{U n ){V4> n + V$ n )e n+1 , V?f)| 

< \ (a(U n )e n+1 V^ n , Vrf)| 

+ |(e n+1 V • {a{U n )V$ n ) + <r([/ n )V$ n • Ve n+1 , 7? n ) | 

< C(|| e n+1 || i2 ||V7f || i2 + ||V^|| L ^||Ve n+1 || L2 || ? f || L2 + llV^ll^llVe^ll^H^llio 



Applying the maximum principle to the elliptic equation (|2.9p shows that ||z,°° < C and 
therefore, 

<c, 

for n = 0, 1, 2, • • • . It follows that 

A(i||e" +1 ||| 2 )+i||Ve^||| 2 

< C\\e n f L2 + C\\e n+1 \\ 2 L2 + C\\q n f m + Cp" +1 ||i 2 

< C\\e n f L2 + C\\e n+1 f L2 + C\\R r l +l \\ 2 L2 , 

where we have used (|3.10p in the last step. By applying Gronwall's inequality, combined 
with (13. 6p . we derive that there exists a small positive constant tq such that when r < tq, 

N 

+ ||^f Hl + £ ll^l^r < Cr 2 . (3.11) 

~ ~ ~ ~ 71=1 

In particular, the above estimate implies that 

\\U n \\ 2 H1 <C (3.12) 

and 

\\D t U n+1 \\ L 2 < ||A^ +1 || i2 + \\D t e n+1 \\ L 2 < C. 
With the above inequalities, we derive from (|2.8p that 

\\U n+1 \\ H 2 < C + C\\V<S> n \\ 2 Li . (3.13) 

Since H 2 (n) ^ C a (Tl) in R d with d = 2,3, ||e n || c - < C. By applying the VF 1 ' 4 estimate 
[61 [21] to U33J), we get 

\\Vr, n \\L*<\Mu n )-v(U n ))V$ n \\ L * 

<C ||e-|| L oo(||vn| L 4 + ||V^|| L 4) 

where Cq is some positive constant. By assuming that Co||e n ||L°o < 1/2, we derive that 

\\Vv n h*<c 
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and (|3.13p implies that 

\\e n+1 \\m < \W n+1 \\m + \\U n+1 \\ H 2 < C (3.14) 

and 

lle^HL- < l|e n+1 ||^ 2 ||e n+1 ||^ 2 < Ct 1 ' 4 . 

From the above derivation, one can see that there exists tq > such that if r < To, then 
Co 1 1 e " 1 1 z, 00 < 1/2 implies Co||e n+1 ||£oo < 1/2 as well as (|3.14p . In addition, we see that 
||V$ n ||i4 < C and therefore, 

max \\U n \\c« < C. (3.15) 

l<n<7V 

By applying Schauder's estimates ([9], page 74) to (|2.9p . we derive that 

max ||V$ n ||c<* < C, (3.16) 

l<n<7V 

which together with (|3.12|) and (|2.9p leads to 

max \W % \\ H 2 < C. (3.17) 

l<n<Af 

Multiplying <^7} by -Ae n+1 and summing up the equations for n = 0, 1, • • • , N — 1, 
we obtain 

N 

max ||e n ||^i + Vr||Ae n ||r 2 
Kn<N n ^— ' L 
~ ~ n=l 

N-l 

< J2 r (lK (J ( C/n ) - ^(^))|V^| 2 ||| 2 + \\a{U n ){V(f) n + V$ n ) • VrfWh + H^Oll*) < Ct 2 . 
n=0 

It follows that 

max |je n || ui < Ct, 

l<n<N 

and 

N N 



^r||AAe n ||| 2 < CT- 2 ^T\\Ae n \\ 2 L 2 <C. 



n=l n=l 

By the theory of elliptic equations [9MH], ll-^* en |li? 2 — C||A-Dte n ||£2 f° r n = !,•■• ,N, and 



so 

A? 



^T\\D t e n \\ 2 H 2 < C. (3.18) 



71=1 



The proof of Theorem 13.11 is complete. I 
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3.2 The fully-discrete finite element solution 

Here we study the error (e^, rfo) of the Galerkin finite element method for the time-discrete 
system (iZ8D- (l2J0l) . 

Theorem 3.2 Suppose that the system (|1.2p - (|1.4p has a solution (u,(p) satisfying (|2.5p . 
Then there exist positive constants ho and tq such that when h < ho and t < tq, the 
fully-discrete system (|2.3|) - (|2.4p admits a unique solution (UJ^,^) such that 

Kh* + Hh* <Ch\ (3.19) 
K\\ H x + H\U<Ch. (3.20) 

Note that the condition of r < to is to ensure that Theorem 13.11 holds. For the given 
Ufi , the error estimate for the equation (12. 4p is given in the following Lemma. 

Lemma 3.1 Suppose that the system I l.ty )- (L4\ ) has a unique solution (u,(j>) satisfying 
fO) . Then 

\\V($ n h -<f> n )\\ L i<C(h+\\e n h \\ L z), 

K - * n \\v> < c{h 2 + \\e n h \\ L2 + ^ 6 IK||i 2 ), 

where (Ufi, <£^) and (U n , $ n ) are i/ie solution of the finite element system 112. 3\) - (2~4\ ) and 
the time-discrete system $2.8\ )- [2~ll^) . respectively. 

Remark 3.1 The proof of the above lemma is similar as that of Lemma 3.2 in [11], in 
which the factor h~ d l® appears when ||eJJ||x,3 reduces to ||e?||_^2 via the inverse inequality. 
More important is that in [11], e^J is the difference between the exact solution of the system 
(|1.2p -( fL~3j) and the fully discrete finite element solution. The restriction for the time-step 
size, r was required when the preliminary error bound ||e^||^2 < C(r + h?) was 

used by induction in the second inequality of Lemma 13.11 However, in our approach, e 1 ^ 
is the difference between the solution of the time-discrete system (|2.8p - (|2.10p and the fully 
discrete finite element solution. Thus, the induction assumption shows that ||e^||^2 < Ch 2 
and then, we can prove the optimal error bound of the scheme unconditionally. 

Proof of Theorem \3.2\ At each time step of the scheme, one only needs to solve two 
uncoupled linear discrete systems. Due to the assumption (jl.5p . it is easy to see that 
coefficient matrices in both systems are symmetric and positive definite. Existence and 
uniqueness of the Galerkin finite element solution follows immediately. It is seen that the 
inequality (|3.20p follows from f|3. 19|) via the inverse inequality. Therefore, it suffices to prove 
(EOT)]) . 

The weak formulation of the time-discrete system (|2.8p - (j2.10p is 

{D t U n+ \ Q + (W n+1 , V£ u ) = (a(C/")|Vd> n | 2 , £ u ), (3.21) 
{a(U n )V^> n , = 0, (3.22) 

for any £^ £ Vh- From the above equations and the finite element system (I2.3p -( [2~4"T) . 
we find that the error function (&h,TjV:) satisfies 

(Aer\ Q + (Ver 1 , V£ u ) 

= {D t (U n+1 - R h U n+1 ), Q + ((a(UJt) - a(U n ))\V<S> n \ 2 , £ u ) 
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+ 2(( ff (E#) - o(U n ))V<S> n • V(*g - $ n ), £„) 

+ (a(^)|V«-^)| 2 ,e«) 

+ 2(a(C/ n )V$ n - V($£-$ n ), e tt ) 

: = (p?+\ u) + (^2 +1 , e«) + e«) + e«) + (3.23) 

and 

{a(UJt)V V l V£^) = -a(C/"))V^, ve ) 

+ (a(E/?)V($" - P£<& n ), V^) (3.24) 

for all Cu,C<p G V h . 

Since r/JJ = on <90, we can take ^ = r/^ in (|3.24p to get 

\\Vrft\\v<C\\<%\\v+Ch, (3.25) 

where we have noted the fact that ||V((/> n — P^(p n )\\L 2 < Ch. With the above inequality, 
from Lemma 13.11 we derive that 

M\\ L 2<Ch 2 + C\\e n h \\ L ,. (3.26) 
Taking £ u = e^ +1 in (|3.23p . the right-hand side is estimated by 

(K? + \ e n h +1 ) < \\e n h +1 \\ 2 L2 + C\\D t U n+l - R h D t U n+1 \\ 2 L2 

< ll^llia + C\\D t U n+l \\ 2 H2 h\ (3.27) 
(P 2 n+1 , e£ +1 ) < C|| e ™ +1 || L2 (||e^|| i2 + \\U n - R h U n \\ L *) 

<C(||e™ +1 ||i 2 + ||e^||| 2 +/ i 4 ), (3.28) 
{RT\ O < C\\el +l \\ L ,{\\e n h \\ L2 + \\U n - R h U n \\ L2 )(\\Vr,l +x \\ L3 + ||*» - P^|| L s) 

< C|| e ^ +1 |Ux(|| e ^|| i2 + C^CHV^Hl. + Cfc) 

< 4e n h +1 f H i + Ce-\\\el\\ L * + ^(CTT^V^^ + Ch) 2 , (3.29) 

(pr 1 , o < cm - * n \\h + cik +i ii| 2 

<C|K +1 ||| 2 +C(/»*+||eJJ||i.) (3.30) 

and 

(i?r\ O < CHeJilli-CllV^lli, + ||V(*» - P^ n )||| 2 ) 
<C/ l " 1 / 2 ||e^ 1 || Hl (K||| 2 +/ l 2 ) 

< e IK +1 H^ + Ce- 1 / i - 1 ||e^||l 2 + Ce^/i 3 . (3.31) 
With the above estimates, by choosing a e small enough, ()3.23|) reduces to 
Adl^llW + llVe^lli, 

<C{\\el +1 \\ 2 L2 + \\el\\\ 2 +h- 1 \\e n h \\h) + Ch 3 + C\\D t U n+1 \\ 2 H2 h\ (3.32) 

which holds for < n < N — 1. 
Now we prove that 

\K\\ L 2 < h 1 ' 2 for 0<n<iV (3.33) 
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by using mathematical induction. Clearly, this inequality holds for n = 0. If we assume 
that this inequality holds for < n < k, then the inequality f|3.32j) reduces to 

A (IK +1 ||i 2 ) + HVOll, <C (\K +1 \\l + \\e n h \\h) + Ch 3 + C\\D t U n+l \\ 2 m h 4 (3.34) 

for < n < k. By applying Gronwalhs inequality, we derive that 

Il4 +1 ||| 2 < Ci/i 3 < h if h < l/C\ /2 . (3.35) 

This completes the induction. 

With (|3.33p . we can apply Gronwall's inequality to (|3.32p and get 

max \\el\\ L 2 < Ch z l 2 (3.36) 

l<n<N 

Since rj? +1 € Hq(£1), from the estimates (|3.25j) - ()3.26|) we see that 

max \\ V n h \\ m <Ch. (3.37) 

l<n<N 

Finally, by applying the W 1,p estimate to the elliptic equation (|3.24|) . we get 

HVr7£ll^ 12 / 5 < - ^IU 12/S H- — ^*-)||^ a:2/s < (3.38) 

and therefore, we obtain a refined estimate: 

(RT\ e n h +1 ) < C||e£|| i6 (||V<|£ 12/s + ||V(* n - P^ n )\\W) 

< e||eg||| 6 + Ce-^V^H^ + Ce" 1 ^^" - P£<^)||£ 12/5 

< e\\e%\\ 2 Hl + Ce~ 1 h 4 . (3.39) 

With the estimates (j3T2T^ - (f3T3TJ| ) and flES]), the equation (pT25D reduces to 

A (K +1 ||| 2 ) + l|V e - +1 ||| 2 < C (||e™ +1 ||| 2 + K||| 2 ) + Ch 4 + C\\D t U n+1 f H2 h\ 
By applying Gronwall's inequality, we get 

max \\el\\ 2 L 2 < Ch 4 . (3.40) 

l<n<N 

The L 2 error estimate of r/^ follows from (]3.26|) and (|3.40p . The proof of Theorem 13.21 is 
complete. I 

Theorem 12.11 follows immediately from Theorem 13.11 and Theorem 



4 Conclusions 

We have presented an approach to obtain optimal error estimates and unconditional stability 
of linearized (semi) implicit schemes with a Galerkin finite element method for the three- 
dimensional nonlinear Joule heating equations. The analysis is based on a new splitting of 
the error into the time direction and the spatial direction, by which the numerical solution 
(or its error) in a strong norm can be bounded by induction assumption and the inverse 
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inequalities without any restrictions on the time-step size. In most existing approaches, a 
time-step condition has to be enforced to bound the numerical solution in a stronger norm. 

Clearly, our analysis can be extended to many other nonlinear parabolic systems, while 
we only focus on the electric heating model in the present paper. A simple example is the 
Joule heating equation with a stronger nonlinear electric conductivity 

a = a(u, Vit) . 

With stronger regularity assumption, optimal I? and H l error estimates in Theorem 12.11 
may be proved without any restrictions on the time-step size r. 

In this paper, we only considered a linear Galerkin finite element approximation. The 
extension to high-order Galerkin finite element methods can be done similarly. For sim- 
plicity, we have assumed that the function g is defined in the domain f2 instead of on the 
boundary d£l. If the function g is defined only on the boundary <9f2, a similar analysis can 
be given by taking the boundary terms into consideration, see |11] for reference. Optimal 
error estimates still can be proved without any condition on the time-step size. 
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